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Abstract 

Negative specific heat is a dramatic phenomenon where processes decrease in temperature when 
adding energy. It has been observed in gravo-thermal collapse of globular clusters. We now report 
finding this phenomenon in bundles of nearly parallel, periodic, single-sign generalized vortex 
filaments in the electron magnetohydrodynamic (EMH) model for the unbounded plane under 
strong magnetic confinement. We derive the specific heat using a steepest descent method and a 
mean field property. Our derivations show that as temperature increases, the overall size of the 
system increases exponentially and the energy drops. The implication of negative specific heat 
is a runaway reaction, resulting in a collapsing inner core surrounded by an expanding halo of 
filaments. 
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I. INTRODUCTION 



While [13| has proven that systems that are not isolated from the environment must have 



positive specific heat, the specific heat in isolated systems can be negative 



Negative 



specific heat is an unusual phenomenon first discovered in 1968 in microcanonical (isolated 



system) statistical equilibrium models of gravo-thermal collapse in globular clusters |12 |. 
In gravo-thermal collapse, a disordered system of stars in isolation under-goes a process of 
core collapse with the following steps: (1) faster stars are lost to an outer halo where they 
slow down, (2) the loss of potential (gravitational) energy causes the core of stars to collapse 
inward some small amount, (3) the resulting collapse causes the stars in the core to speed up. 
If one considered the "temperature" of the cluster to be the average speed of the stars, this 
process has negative specific heat because a loss of energy results in an increase in overall 
temperature. 

In the intervening four decades, negative specific heat has been observed in few other 
places. In a magnetic fusion system or other thermally isolated plasma, should negative 
specific heat exist, the related runaway collapse could have profound implications for fusion 
where extreme confinement is critical to a sustained reaction. 

Our results have general applicability to vortex systems. However, in this paper, we 
address a plasma model known as the electron magnetohydrodynamical (EMH) model, where 
we report finding negative specific heat. Our findings are based on a mean- field approach 
to the statistical equilibrium of the system. 

Typically, magnetohydrodynamic plasma models are two fiuid models, requiring equa- 
;ions governing the electron motion and equations governing the ion motion coupled together 
1^. The EMH model bypasses the two-fiuid model by representing the electron fiuid and 
the magnetic field as a single, generalized fiuid with a neutralizing ion background that is 
stationary on the timescale chosen. 

The EMH model takes the magnetic field, B = V x A, and the charged fiuid vorticity, 
a; = V X V, and combines them into a general vorticity field i7 = V x p where the generalized 
momentum, p = mv — eA, m is the electron mass, — e is the electron charge, v is the fiuid 
velocity field, and A is the magnetic vector potential field. For a brief overview of the model, 



see 



16| . A detailed model discussion can be found in 6| 



Our goal is to find an explicit, closed-formformula for the specific heat of this vortex model 
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in statistical equilibrium given an appropriate definition for energy and a microcanonical 
(isolated) probability distribution. Our approach is to describe the statistical behavior of a 
large number of discrete, interacting vortex structures and consider the limiting case. We 
hypothesize that the specific heat is negative. 



II. QUASI-2D VORTEX MODEL IN STATISTICAL EQUILIBRIUM 

In this section we define nearly parallel vortex filaments which are the "particles" of our 
model. They are asymptotically thin vortex tubes that are polarized to be nearly parallel 
to the axis of rotation (or magnetic confinement; in the case of the EMH model, both can 
be considered to be present). We go on to give (but not derive) the equations of motion 
and Hamiltonian for these filaments and the statistical distribution for the filaments in a 
microcanonical distribution. 

First we define nearly parallel vortex filaments: 

Definition II. 1 Nearly parallel vortex filaments are smooth curves with a complex param- 
eterization tpilr^t) where ipi{T,t) = Xi{T,t) + iyi{T,t) and r G [0, 1], and t is time. They are 
periodic, ipi{Q,t) = ipi{l,t). 

If we take any two values of t, tq and ti such that tq < ti, and let At = ti — tq such 
that At G 0(e) where e ^ 1, then for any filament i, the amplitude is of order , i.e. 
mri)-Ur^)\eO{e^). 

In words this means that, for a small rise of length e in the filament, the amplitude must 
be on the order of e^. This assumption guarantees a certain degree of straightness in the 
filament that allows for the derivation of the quasi-2D equations of motion. 

The N coupled, non-linear Schrodinger equations for the motion of these curves are the 



following 



dt Ot^ — Tp 

The PDE leads to a convenient energy functional 

iV -, o , / \ 2 ^ N N 



^ J2J2^og\U^)-^P,iT)\dT, (2) 

" i=l i=l 7=1 
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where a is the core-structure constant in units of energy /length. The first term in the energy 
derives from a local-induction approximation (LIA) that causes Brownian variations along 
the length of filament in the plane. The second term is the coupling term and results in 
repulsion between filaments, a typical 2D Coulomb interaction that happens only within 
each plane and not between planes. 

The energy is the first conserved quantity. The second conserved quantity is angular 
momentum, 

N 1 

i=l 

We define the enthalpy to be — +pM]sr, where p is a fixed parameter that we call 
pressure. Because energy and angular momentum are conserved, enthalpy is conserved. 

We define the micro-canonical probability for the statistical equilibrium of the system of 
N curves. First 

Definition II. 2 Let the set of states U be the space of all sets of N smooth, complex func- 
tions on the interval [0, 1], u = {'4'i}i=i...N, such that J \ipi{T)\'^dT < oo. 

The functions in U include both nearly parallel and not-near ly parallel functions. How- 
ever, we will argue that not-nearly parallel functions have negligible effect on the statistics. 

Definition II. 3 Let : C/ — > [0, oo) he a Gibbs density functional and Pjv : U ^ [0, 1] be 
a Gibbs probability density defined by the equations 

Gn(u) = 5(NHo - Hn)5(NR^ - Mn) (4) 

„ / X Ga[{u) , , 

Pn{u) = (5) 

where = j GN{u)du is a normalizing factor called the partition function, and S is a 
Dirac-delta function. 

The constants Hq and are the enthalpy and angular momentum per vortex filament per 

period respectively. 

The density and probability density functionals define the intersection of the function- 
space "areas" of the enthalpy-surface and the angular-momentum-surface. If the enthalpy 
surface is the level set Hsurf — {u\Hn — NHq} and the angular momentum surface is 
Msurf — {u\NR'^ — Mn}, then the intersection of the two sets is A — HgurfClMsurf: and 



Zn is the size of the set A. Any probabihty p = PN{u)du describes the size of the 
intersection of B and A normahzed by the size of A. The probabihty density describes the 
statistical equihbrium of the model. 

Note: An equivalent distribution, the intersection of the energy surface and the angular 
momentum surface, is less appropriate in this case because it does not make the pressure 
constant, p, explicit. Also note that the angular momentum is an output parameter not 
an input. Again, we include it to make it explicit. 



III. THE MEAN-FIELD APPROXIMATION 



To find an explicit, closed- form formula for the specific heat, we need to have an explicit, 
closed- form formula for the partition function Z^. Current mathematical understanding 
makes a direct integration impossible, but an approximation can remove the difficulty. This 
approximation is called a mean-field approximation because it averages the effect of all 
the filaments on one and replaces the other filaments with an external field based on the 
average. Because the interaction term is based on the distance between each of the filaments, 
we replace, at some unknown cost, this distance with an average value. In the next section. 
Section IIII At we discuss the cost of this assumption. 

Theorem III.l Assuming that for any given filaments i and j and plane r the random 
variables ipii'^) ^'^c uniformly distributed on a circle of radius 2R, then 

N N „i 



E 



where E denotes expectation value. 



— log/^^ 



(6) 



The expectation value can be brought into the sums and integrals with no problem, 

N N '] ^ N N 

4 



E 



N N I 1 N N I 

4 EE / logl^.-^.Prfr =4EE/ E[log\^,-^P,\']dr. 

i = l j = l "'0 J j^;^ Jo 



(7) 



Now let ipilr) = zi and ipjir) = Z2, and consider the expectation in polar coordinates, 

'•271 r^TT r2R r'iR 



E [log \zi - 2;2p] 



{-kAR 



2^2 



^0 ^0 ^0 



log \r'l + rl — 2005(6*2 — ^i)] ridrir2dr2dOid02. 



(8) 
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This integral has been shown to evaluate to \ogR^ plus a constant that we can drop 
without loss of generality . 

This concludes the mean-field approximation. The energy functional now reads: 



E' 



N 



a 



1 N 

El 



i=l 



dr 



dr 



— logi?'. 



(9) 



Each of our previously defined functionals, Hjy, Gjy, Pjy, and Z^, now has a mean-field 
version, and, for the sake of simplicity of notation, we drop primes and refer only to the 
mean-field functionals. 



A. A circle of radius 2R 

The assumption of uniformity is not as drastic as it may appear. First of all, just because 
we assume a uniform distribution for the interaction energy does not make the resulting 
distribution uniform. In fact, the distribution is only uniform if the self-energy is zero. 
Therefore, the assumption is not that Pn is uniform but that assuming that it is uniform 
for the sake of gaining a simpler interaction term does not change the physics significantly. 
There is significant justification for this assumption. 

We argue using results from related statistics. In a previous paper using Monte 
Carlo simulations, we showed this same mean-field approximation to be extremely effective 
in the canonical (non-isolated) case of nearly parallel vortex filaments with distribution 
= exp{—i3H]\r)/Z'f^, where = exp(— /3iJjv) and /5 is the inverse temperature 
parameter. The canonical distribution and the micro-canonical distribution are often equiv- 
alent, and, even when they are not, their distributions can have many of the same properties. 
As we will show, with the mean-field approximation in place, the formula we obtain for i?^ 
in terms of the parameters from the micro-canonical distribution, P/v, is identical to the 
formula we obtained in our previous paper for the canonical distribution, P^. 

NB: If the distribution is uniform with radius 2R, then the mean angular momentum 
must be NR^, and, if the distribution is uniform and the angular momentum is NR^, then 
the radius of the distribution must be 2R. This is from the definition of M^- 
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IV. MAXIMAL ENTROPY 



The specific heat formula that we are about to derive is the specific heat of the most- 
probable macrostate. A macrostate in our case is specified by the number of particles, 
A^, the enthalpy, Hq, and the pressure, p and is a collection of microstates: A{N, Hq, p) C U. 
There are many macrostates for each fixed set of parameters, and each has a different total 
entropy, S, temperature, T, angular momentum, i?^, specific heat, Cp, etc. However, only 
one of these macrostates is likely to occur, the one with the largest entropy. To see why 
note the formula for entropy-density by Shannon [l^ . 

SNiu) = -PNiu)\ogPNiu), (10) 

where u E U is a microstate. The bigger the probability, the bigger the entropy. Therefore, 
the macrostate with maximum entropy is also the most-probable: 

Definition IV. 1 The most-probable macrostate is a macrostate Ump C U such that the 
entropy of Ump is maximal. The maximal entropy is defined as 3^°"^ = J Si\f{u)du, where 
the integral is over u G Ump- 

This definition is extremely valuable because it gives us a way of describing state variables 
such as temperature and specific heat with fixed values, rather than the true, fluctuating 
ones: 

Definition IV. 2 The inverse temperature of the system is defined to be the change in 
maximal entropy with respect to total enthalpy, 

1 Q^max 

^' = T = d{NH^y ^^^^ 

and the specific heat ( at fixed pressure p ) the change in total enthalpy with respect to tem- 
perature at maximal entropy, 

djNHo) ., d{NH,) 



Defining these variables in terms of the now fixed maximal entropy (states of greatest dis- 
order) and fixed enthalpy makes temperature and specific heat fixed for a given set of 
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parameters. Redefining state variables that are not fixed in reality into those that are is the 
essence of statistical mechanics. 

With these definitions, we can begin to calculate the formula for the maximal entropy 
which will lead to specific heat in the next section. We now limit our investigation to only 
the most-probable macrostate: 

Theorem IV. 3 Assuming that all states are in the most-probable macrostate, i.e. Ump = 
U , in the limit as N oo and with the necessary, non-extensive scalings (3q = (3qN , 
H'q = Hq/N, a' = a/N, and p' = p/N , the maximal entropy per filament is 

^-« + T">^(«^)-25^-''i'''«- (13) 

where 

^ = 8^' • ^^^^ 

We proceed to give a brief proof of the formula: 

We can define the maximal entropy in terms of the partition function: S^"-^ = logZ^v 
using Equation [10] above. We re-write it to look like this: 

e^^ = J Dij5iNHo- Hn)S{NR^ - Mn), (15) 

where NS = S'^"'^ and S is the maximal entropy per filament, and / Dip is a Feynman 
integral operator 5|]. When we take in infinity, we will be left with S rather than S^°-^. 

In order to apply steepest-descent methods, we replace the Dirac-delta function with its 
integral representation. 



f3o+ioc 

27Ti' 



where jSo is defined to be the /3-value at which the integrand attains its maximum value. 

The integral is clearly finite because the exponential is quadratic negative definite, and 
the remaining delta function is only a constraint making the domain of integration smaller. 
One could argue that Fubini's theorem does not apply to functional integration, but, in the 
case of this particular, traditional Feynman integral, it does apply, and one could easily 
show it using Feynman's piecewise linear segment approximation and taking the limit as the 
number of segments becomes large jsl. The order of integration can be rearranged to give. 



/ ^e^^^« / D#(iVi?2_M^)e-^^-. (17) 



This rearrangement has produced an expression that is famihar from our previous papers 
I2] and fl, namely the canonical partition function, Zcan = J D%l)5{NR'^ — Mm)g l^^^ . 

In those papers, we showed that if F = — limTv^oo ;g7v -^^can) given the scaling that 
(3' = pN, a' = a/N, and p' = p/N, then 



F = p'R'^ -l/ilogR^ + ^ 



where 

2 _ f^'^a' + ^/p"^a'^ + 32a'p'^p' 
~ 8a'/3'V ' 

In order to use our previous results, we need to take the limit (and scalings) on Equation 

S= hm llog T"^'" ^e^'^^oz^. 

The steepest-descent argument from and |7| says we can replace all the instances of 
with /3o and lose the integral over 



S = 13'qH'q + limAT^oo log 



which proves the formula. 



V. AN ASTRONOMICAL ANOMALY: NEGATIVE SPECIFIC HEAT 

Knowing the maximal entropy for any system for a given set of parameters, allows one 
to determine any of its state variables. However, in our case we have derived a maximal 
entropy that is dependent on a state variable, temperature. Therefore, we cannot determine 
the specific heat until we can define the temperature in terms of the input parameters. 

We find the unknown multiplier, /3o, by relating the enthalpy per filament parameter. Ho, 
to the mean enthalpy: NHq = {Hn) [TI]. 

Theorem V.l The enthalpy, H'q, for the infinite-N system (as defined above) is related to 
the inverse temperature for that system by the following: 

^ f Po 1 d2 , -'- I a/ / r>2 
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By definition the average enthalpy is given by 

, , ^ / DijHr,6{NHo - Hn)6{NR^ - M^) 

^ J D^6{NHo - Hn)6{NR^ - Mm) ' ^ 

which has the following integral representation: 



(Hn) 



'/No- 
where we have replaced Hj~^ exp(— /5i/jv) with —d/d(3 exp(— /S/Z^v). Now in order to continue 
we need to bring the integral over Dip inside the derivative so that we can use our previous 
steepest-descent results to simphfy the expression. We prove that this is allowed in Appendix 

m 

Let H = N~'^Hn- Applying the limit and scalings as before (dropping the limits of 
integration), we exchange the derivative and functional integral: 

1 / 'S^'''"'^ (-W I Di,e-^'^-5{NR' - M^] 

{H) = lim 



jv^oo N J f^e^'^^o / D^e-P'"^6{NR^ - M^) 

Again according to the steepest-descent argument in 71], only values of the integrand at 
/3q have any contribution to the integral as N becomes large: 

1 (- / D^e-^'o^'^diNR^ - Mn) 
~ iv / Dilje-P'o^NS^NR^ - Mn) 
The fraction simplifies to give 

{H) = lim^^oo (- log / D^e-f'oH^diNR^ - M^i 

= w/oF. (20) 

where F is defined in Equation [TSl (Note that direct evaluation can show that switching 
the limit and derivative is admissible.) Since H = Hq, the formula is proven. 

Based on the definition in Def. IIV.2I and the two theorems, i.e. the formula for maximal 
entropy and the formula for mean energy and consequently temperature, a Corollary is the 
formula for specific heat. 

Corollary V.2 The specific heat has the form: 

c,^^4( ^ -iV (21) 
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A little algebra is all that is required. Therefore, like its astronomical sibling, this vortex 
system has negative specific heat, a hallmark of meta-stable states (see next section for 
details). We would also like to point out that the key to obtaining negative specific heat is 
the non-extensive nature of the system. Extensive systems, by mathematical proof, cannot 
have negative specific heat pL5,]. Therefore, the scaling that we introduce to make the limit 
non-extensive is what breaks this proof's assumptions and allows for negative specific heat. 

VI. RUNAWAY REACTION: A THOUGHT EXPERIMENT 

All of these mathematics rely on the mean-field assumption, which we argued was reason- 
able in Section IIII A[ We cannot tell if the specific heat of the original mathematical model 
is negative directly, but consider that all we have done is taken a complicated expression 
for the interaction and reduced it to something that is intuitive, namely that as the system 
expands in average size the interaction increases logarithmically. We offer the following 
thought experiment to demonstrate how negative specific heat causes the system to behave: 

Suppose that we have the system at an enthalpy Hq and decrease the enthalpy to Hi < 
Hq. Several possible corrections can occur in each of the following enthalpy terms: self- 
energy, interaction energy, and angular momentum. Either the filaments become straighter, 
decreasing the self-energy, or they can move apart, decreasing interaction, or they can move 
closer, decreasing angular momentum. Most likely, the correction will be a combination of 
the three depending on which maximizes the entropy. There are two sources of entropy: (1) 
increased misalignment of the filaments (decreased straightness) (2) expansion of the system 
as a whole. Some balance of the two will occur. 

Now suppose the correction proceeds as follows: the filaments closest to the origin squeeze 
together. They become straighter and decrease the self-energy, increase the interaction 
energy, and decrease angular momentum. The filaments further away from the origin move 
outwards decreasing the interaction energy but increasing the angular momentum and the 
entropy related to expansion. They also increase the entropy by becoming less straight, 
increasing the self-energy. The three sources of enthalpy are balanced so that the total 
enthalpy is reset to Hi due to the decreased angular momentum and self-energy in the 
origin and the decreased interaction energy in the outer halo. 

The overall entropy should decrease with decreased enthalpy because we have positive 
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temperature. The system expanded, but the straightness decreased at the origin, so let 
us say that the total entropy decreased. (Of course, in a point vortex system, this cannot 
happen, but here it is possible.) Because pressure p is constant and the volume in the 
center decreased, the temperature there increased. Meanwhile, the filaments that moved 
outward saw an increase in volume and corresponding decrease in temperature. Because the 
cost of moving outward in terms of increased angular momentum is with the square of the 
distance, the ones moving outwards will not move outwards as much as the ones moving 
inwards. Therefore, the overall temperature increases with the decreased enthalpy. This 
same effect is observed in globular-clusters Q,H. We propose that this effect will show 
up in a numerical simulation of the micro-canonical system for certain parameter regimes. 



VII. CONCLUSION 



We have shown a way of simplifying the interaction of a nearly parallel vortex filament 
model with an intuitive mean-field approximation and calculated a formula for the specific 
heat at maximal entropy for the system in isolation. The specific heat was shown to be 
negative indicating that the vortex system is meta-stable. With a thought experiment we 
proposed what the negative specific heat would mean, i.e. that vortices would separate out 
into a core and halo and that the core would fall in on itself. This core would likely result 
in vortex merger creating a large "hole" at the origin with a large vortex swirl around it 
and many smaller filaments surrounding it. Although this work is entirely theoretical, we 
propose that its results can and will be observed in computational and experimental settings. 

In the introduction we suggested the EMH model as an application because of its clear 
relationship to magnetic nuclear fusion in which confinement is the key to attaining a self- 
sustaining fusion reaction. The runaway collapse we propose may lead to such a reaction. 



APPENDIX A: PROOF OF INTERCHANGE OF DERIVATIVE AND FUNC- 
TIONAL INTEGRAL 

The following lemma relies on what is known as the broken segment or piecewise linear 

n n 

approximation of the Feynman paths ([5|,[10|). Because functional integrals require different 
and less familiar mathematical machinery than ordinary iterated integrals, the broken seg- 
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ment model provides a way to prove things about functional integrals of our particular type 
(rather than general functional integrals) using ordinary and familiar integral theorems (e.g. 
Fubini's). To get the broken segments, we take the complex function ipi'^) and approximate 
it with a piecewise linear vector, \E', of length M, i.e. \1' = {ip{Ti),tp{T2), . . . ,iP{tm)) where 
Tj+i = Ti + l/M and tm+i = Ti- The functional integral operator f Dip = limjv/^oo / d'^/a^ , 
where a is a necessary scaling factor. 

Lemma A.l The following equation is true: 

lim r d<i'/a^^S,e-'^^^^^'^'>5(NR^ - Mj^(M)) 

M^oo ' op 

= ^exp(-/5ArF), (Al) 

where a = 7i / (aPoM) . 

There are two separate issues here: 

1. whether the derivative can be brought out of the integral: 

jd^/a^ ^e-^^^(*^)5(Ari?2 _ M,v(M)) 

= ^ / ci^/a^e-^^^(*^(5(Ari?2 _ Mjv(M)), (A2) 

2. whether the derivative can be brought out of the limit: 

Jim ^ / rf^/a*^e-'3^^(^'^)5(iVi?2 _ Mjv(M)) 

= ^ limM^oo / rf^/a*^e-'3^^(^^)(5(Ar/?2 _ Mn{M)). (A3) 

Equation IA2I is true if and only if 

c/^/a*^e-^^^(*^)5(A^i?2 _ MNiM)) < oo, (A4) 



and the integrand is differentiable. The equation is integrable because the integrand is 
positive definite (and positive). Therefore, it is a simple Gaussian integral. The smoothness 
of the integrand guarantees differentiability. 

For IA31 we prove it explicitly, i.e. calculate the derivative first and then take the limit. 
If M is large, we can make the approximation: 



/ 



(i^/a^e-^^^(^)5(Ari?2 _ Mn{M)) ^ e'^^^^, (A5) 
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where 



and 



F(M) =pB? - NH\ogB? - !^RHvo - ^) - J l»8('!o + « " 1)') (A6) 



Taking the derivative, 

d 



^^g-/3iVF(M) ^ NF{M)e-^''^^^^Ff,(M), (A8) 



it is trivial to show that 



hm 7VF(M)e-^^-^WF^(M) = NFe-^''^Fp = -^e"'^^^. (A9) 
M— >oo a/3 
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